Filamentation instability in a quantum plasma 
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The growth rate of the filamentation instability triggered when a diluted cold electron beam passes 
through a cold plasma is evaluated using the quantum hydrodynamic equations. Compared with a 
cold fluid model, quantum effects reduce both the unstable wave vector domain and the maximum 
growth rate. Stabilization of large wave vector modes is always achieved, but significant reduction 
of the maximum growth rate depends on a dimensionless parameter that is provided. Although 
calculations are extended to the relativistic regime, they are mostly relevant to the non-relativistic 
, one. 
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I. INTRODUCTION 



There is currently a growing interest in quantum plasmas in connection with many areas of physics like small 
electronic devices, astrophysics or laser plasma interaction (see Refs. [H, HJ and references therein). As far as plasma 
i instabilities are concerned, the quantum theory of the two-stream instability has been elaborated in Refs. [H, |j, [H, 
while quantum effects on the filamentation instability are still to evaluate. The purpose of this work is to deal 
with this issue which could be relevant for some astrophysical settings, or the Fast Ignition Scenario for Inertial 
Confinement Fusion [7j where a relativistic electron beam is supposed to reach the partially degenerate dense core of 
03 , a pre-compressed target. 

We thus consider an infinite and homogenous cold relativistic electron beam of velocity Vj,, relativistic factor 7& 
and density rib entering a cold plasma of electronic density n p . Ions form a fixed neutralizing background of density 
Tib + Tip and the beam prompts a return current in the plasma with velocity V p such as n p "V p = ribVb- In a first 
attempt to quantity quantum effects for the filamentation instability, we implement the quantum fluid formalism 
relying on the conservation equations for the beam (j — b) and the plasma (J = p) , 

^ + V-(n iVi ) = (1) 
and the force equation with a Bohm potential term 0], 




> 

?| ^ + (v,-V)p J =- (Z (E + ^-^) + £-V(-^], (2) 

where q > and m are the charge and mass of the electron, nj the density of species j, and pj its momentum. 
Quantum corrections are clearly contained within the Bohm potential, which can be obtained from the moments of 
j the non-relativistic Wigner function [2j. The present calculations are here extended to the relativistic regime, but a 
correct treatment of quantum effects in that limit should rely on the moments of a relativistic Wigner function such 
as the one described in Refs. [1, Nevertheless the modified force equation © has to our knowledge no simple 
relativistic counterpart yet, and present relativistic results make sense physically. This is why we found relevant to 
deal with the relativistic regime with the present equations. Although the following theory is only correct in the 



non-relativistic case, it could also be interesting to compare future relativistic results with the present ones. 
Assuming quantities perturbed according to exp(ik • r — itot), we find the linearized equations, 



_k_Vj! 

to — k • Vjo 



n jl = n j o- — - , (3) 

and 
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where subscript and 1 stand for the equilibrium and perturbed quantities, and jj for the relativistic factor of the 
unperturbed beam and return current. Note that the diluted beam hypothesis <C n p implies a non- relativistic 
return current with j p = 1. 

We now choose to align the beam velocity along the z axis and k along the x axis, and calculate the dispersion 
equation following a quite standard procedure [10|. Expressing the current J in terms of the electric field E and 
inserting it into a combination of Maxwell's equations, 



^kx (kxE)+E + J 



0, 



allows for the derivation of the dielectric tensor. Its determinant finally gives the dispersion equation 



(5) 
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where uj p is the electronic plasmas frequency, and 



e 



with O c = (hujp/2mc 2 ) 2 . Although this quantity remains small even in a very dense plasma with G = 1,3.10~ 7 
assuming n p = 10 26 cm -3 and V b ~ c, we will check in the sequel that quantum effects may not be negligible for the 
filamentation instability. 



II. QUANTUM EFFECTS ANALYSIS 

Let us start plotting on Figure [T] the growth rate in terms of Z for a diluted beam passing through a dense 
plasma, in the non-relativistic and relativistic regimes. We make the comparison between the calculations for C = 
and O c — 1, 3.10 -7 . In both cases, it is straightforward that although the maximum growth rate is not modified (and 
we will see why later), quantum effects introduce a cut off at large Z whereas the growth rate just saturates in the 
classical limit = 0. 

The cut off value of Z can be calculated exactly in terms of O noting that a vanishing growth rate implies that Eq. 
^ be verified with x — since the real part of the root is for the filamentation instability. The resulting expression 
can be solved exactly in terms of Z and gives the cut off wave vector Z m , 



(9) 



Accounting for a < 1, we find 



It turns out that the parameters 4a/Of3 2 measures the strength of quantum effects as far as Z m is concerned. For 
the sake of physical interpretation, let us consider the most common case where this parameter is large with respect 
to 1, and set ~ 1. Equation (TTU|) can be developed and yields 
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FIG. 1: Growth rate 8 in ui p units of the filamentation instability for a diluted non-relativistic (upper graph) and relativistic 
(lower graph) beam. The plasma density in both cases is n p = 10 cm" " 3 yielding 9 C = 1, 3.10" 7 . The range of unstable wave 
vectors is dramatically reduced whereas the most unstable mode and its growth rate are only weakly affected. Parameter Q 
given by Eq. (|14p is here 307 and 439 respectively. 



which is the inverse of the beam quantum wavelength. Simply put, perturbations which scale length is smaller than 
the beam quantum wavelength are relaxed and do not grow. Also remarkable, the 7& dependance is lost between Eqs. 
© and (flT)]) because of the diluted beam hypothesis, so that it can be said that to a large extent, quantum corrections 
to the instability domain in the relativistic regime do not depend on the beam energy as long as they are described 
with the correction term inserted in Eq. ©. 

While quantum effects on Z m are necessarily important since Z m — oo in the classical case, effects on the largest 
growth rate are even more relevant for the physics of the system. If we recall that in the classical limit, the growth rate 
8 saturates for Z > (3 with S ~ /3y/ a/jb the occurrences of in dispersion equation ((6| indicate that quantum 
effects near the most unstable wave vectors are negligible if 

e>z 4 \ z ^p < x 2 \ n r —, 

' ' 'x~f3yf a/71, ' 

QZ*\ z „b < j b x 2 \ . ri -. (12) 
These two conditions are clearly met if the first one is satisfied, 

9/? 4 « P 2 a/ lb & -^-2 » 1. (13) 

We here recover the a/O/3 2 parameter already evidenced when dealing with the largest unstable wave vector, but we 
find it must dominate the beam relativistic factor if quantum effects are to be negligible with respect to the most 
unstable mode. Therefore, depending on the inequality (|13p . we will find situations where quantum effects stabilize 
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FIG. 2: Region of the (a, /3 or jb ) plan where fi ^> 1 for n p — 10 25 and 10 30 cm 3 . The curves picture the border Q, — 1 and 
the shaded areas indicates the side of the border where quantum effects are negligible with S> 1. 
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the instability at large Z while changing, or not, the most unstable mode. The parameter 

Q^fbP 6 C 76 

is thus the most relevant one from the physical point of view because a change in the fastest growing mode brings 
more physical consequences than the stabilization of the non-dominant large Z modes, fl 3> 1 implies few quantum 
effects for the most unstable mode whereas fi 1 means strong quantum effects. 

The physical interpretation of fl can be elaborated after the one we used for Z m in Eqs. (|10lll[) . We found that 
wavelengths larger than the beam quantum wavelength are stabilized. Let us now construct the ratio of k m as given 
by Eq. (jTTJ) with the plasma skin depth uj p /c wave number, where the instability reaches its maximum in the classical 
regime. We take the fourth power of this quantity and consider the relativistic beam regime and find, 

We thus check that once added a relativistic contraction factor 1/jb, the parameter £1 can be constructed from the 
ratio of the plasma skin depth X s to the beam quantum wavelength k^ 1 . When A s is much smaller than fc~ 1 , the 
perturbations which should give rise to the most unstable modes are "blurred" by quantum effects, resulting in a 
tamed instability. 

In the non-relativistic regime, fl = a/3 2 /0 C can be traced back, after some manipulations exploiting the identity 

n b /n p = wl/ul = Vp/V b , to 

(16) 

where A e = h/mVb is the de Broglie length of a beam electron. Quantum effects are here important for the maximum 
growth rate when the beam electron de Broglie length is much larger than the plasma skin depth. 

We picture on Figure [2] the region of the (a,/3 or 75) plan where Q 3> 1 for two different plasma densities. From 
Eq. (T4"|) and Fig. [21 we infer that quantum effects tend to be important for beams with j3 <C 1 or 7b ^> 1, that is, 
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FIG. 3: Growth rate 5 in uj p units of the filamentation instability for a non-relativistic (upper graph) and a relativistic beam 
(lower graph) and Q, <C 1. Here the most unstable mode and its growth rate are significantly modified by quantum effects. 



the non-relativistic and ultra-relativistic limits. The smallest beam densities are required in the intermediate regime 
for 7 b = \/2 with a < 26 c . 

Figure [T] already displays two systems with Q, = 307 and 439 ^> 1. As expected, according to the previous analysis, 
the maximum growth rates are weakly affected by quantum effects whereas the ranges of the instability are reduced. 
Figure [3] pictures now one non-relativistic and one relativistic situations with small f2 parameters. The maximum 
growth rate is significantly reduced in both cases so that quantum effects should have easily observable consequences. 
Noteworthily, the relativistic scenario requires a much higher plasma density than the non-relativistic one, where f2 
can be lowered down through the 1 factor at the numerator. Quantum effects in the relativistic regime should thus 
be expected mostly in astrophysical settings, as already noted in the context of quantum magnetohydrodynamics Q 
for example. 



III. DISCUSSION AND CONCLUSION 



We have implemented the quantum fluid equations to evaluate quantum corrections to the filamentation instability. 
Quantum effects stabilize the instability for large wave vectors and can significantly reduce the most unstable mode 
if parameters fl defined by Eq. (fl"4"l) is much smaller than unity. Under such conditions, quantum effect can reduce 
the filamentation instability. We could check that there is eventually no quantum cancelation of the instability, 
fndeed, a numerical exploration of the regime Q <C I shows that the maximum growth rate behaves like VTi in this 
limit. Given the requirements on the parameters for quantum corrections to become significant, the necessity of a 
quantum treatment may be restricted to some astrophysical scenarios when relativistic beams are involved. Even in 
the conditions of the Fast fgnition Scenario for fnertial Confinement Fusion, cold fluid quantum corrections are only 
about 10%, as n p — 10 26 cm~ 3 , rib = 10 21 cm~ 3 and jb = 5 yields £1 — 159. Nevertheless, a correct relativistic 
treatment should rely on the moments of the relativistic Wigner function, so that present results are only valid in the 
non-relativistic case while their relevance for the relativistic regime may be only qualitative. 

Let us add that it would be very interesting to investigate the joined effects of temperature and quantum degeneracy. 
In the classical relativistic limit, it is well known that kinetic effects can literally suppress filamentation [H|, [l3j]. Since 
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quantum effects can also have some stabilizing function, a quantum kinetic treatment of the problem could well 
unravel even more pronounced stabilizing quantum effects than those investigated here. 

In the same way, it has been known for lon g th at filamentation instability can be significantly reduced, if not 



completely suppressed, in a magnetized plasma 14j . While various quantum effects have been recently investigated 



in the case of a magnetized plasma [HI, [3, [T3| , the behavior of the filamentation instability in such settings should 
be a subject of much interest. 
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